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ABSTRACT 

We report precision Doppler measurements of three intermediate-mass sub- 
giants obtained at Lick and Keck Observatories. All three stars show variability 
in their radial velocities consistent with planet-mass companions in Keplerian or- 
bits. We find a planet with a minimum mass Mp sin i = 2.5 Mj up in a 351.5 day 
orbit around HD 192699, a planet with a minimum mass of 2.0 Mj up in a 341.1 day 
orbit around HD 210702, and a planet with a minimum mass of 0.61 Mj up in a 
297.3 day orbit around HD 175541. Mass estimates from stellar interior models 
indicate that all three stars were formerly A-type, main-sequence dwarfs with 
masses ranging from 1.65 M & to 1.85 M . These three long-period planets would 
not have been detectable during their stars' main-sequence phases due to the 
large rotational velocities and stellar jitter exhibited by early-type dwarfs. There 
are now 9 "retired" (evolved) A-type stars (M* > 1.6 M ) with known planets. 
All 9 planets orbit at distances a > 0.78 AU, which is significantly different than 
the semimajor axis distribution of planets around lower-mass stars. We examine 
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the possibility that the observed lack of short-period planets is due to engulfment 
by their expanding host stars, but we find that this explanation is inadequate 
given the relatively small stellar radii of K giants (i?* < 32 R Q = 0.15 AU) and 
subgiants (i?* < 7 R Q = 0.03 AU). Instead, we conclude that planets around 
intermediate-mass stars reside preferentially beyond ~0.8 AU, which may be a 
reflection of different formation and migration histories of planets around A-type 
stars. 

Subject headings: techniques: radial velocities — planetary systems: formation- 
stars: individual (HD 192699, HD 210702, HD 175541) 



1. Introduction 



Very little is known about the occurrence rate and orbital properties of planets around 
A-type stars, corresponding to stellar masses ranging from 1.6 M & to 3.0 M & . Inspection of 
the Catalog of Nearby Exoplanets (CNE)0 reveals that only 6 of t he 173 stars with s ecurely 
detected planetary companions have masses in excess of 1.6 M jButler et alibood l This 
small number of detections is not a true reflection of the occurrence of planets around A-type 
stars, but rather the result of a strong selection bias against early-type, main-sequence stars 
in precision Doppler surveys. 

Measuring precise Doppler shifts of early-type dwarfs is complicated by their rotation- 
ally broadened spectral fea tures, high surface temperature s, and high level s of ex cess radial 
velocity noise, or "jitter" JSaar et all fl998l : IWrightl I2005T ). iGalland et all J2005h find that 
Doppler precision for early-type dwarfs is limited to ~ 40 m s _1 at spectral type F5V, and 
90-200 m s _1 for A stars, rendering Doppler measurements of these stars sensitive only to 
planets with large masses and short orbital periods. The lowest mass companion so far de- 
tected around an A star is the brown dwarf orbiting HD 180777 ( Galland et al.ll2~006l ). Even 
though the 28 day orbital solution has a large velocity semiamplitude, K = 1200 m s -1 , the 
signal is only a 3<r detection above the stellar jitter and measurement uncertainties. 

Most of what is known about planet formation around intermediate-mass stars comes 
from two primary sources: direct imaging of disks around young stars and Doppler detections 
of planets around evolved stars. While A-type dwarfs are poor Doppler targets, their high 
intrinsic luminosities facilitate the detection and direct imaging of material in their circum- 
stellar environments. More than a decade before the discoveries of the first extrasolar planets, 



For the updated catalog of extrasolar planet and their parameters see 
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evidence of planet formation outside of our Solar System came from the infrared detection o f 
collision-ge nerated dust around th e A-type, main sequence stars Vega (lAumann et al.lll984l ) 
and j3 Pic ( ISmith &: Terrild Il984l ) . Since then, advances in high-contrast imaging have re- 
sulted in the detect ion of an optically thick disk around a pre-main-sequence Herbig Ae star 
( jPerrin et al.ll2006l ). as well as scattered light images of optically thin "debris disks" around 
11 main-sequence stars — the majority of which have spectra l types F5V or earlier (Table 2 of 
Kalas et al.ll2006l ; ISchneider et al.ll2006l ; IWahhaj et al.ll2007l . and references therein). Recent 
observations of the debris disk around the young A star Fomalhaut have revealed a pertur- 
bation in the disk structure that may be due to the influence of an orbiting Jovian planet 
( IKalas et al.l 120051 ) . Studying the relationships between the architectures of disks around 
young A stars and the distribution of planet properties around their older counterparts will 
provide key tests of planet formation models. 

A key to finding planets around A stars using Doppler methods is provided by the 
effects of stellar evolution. As stars evolve away from the main sequence, they become 
cooler a nd rotate slower, wh i ch increases the numb e r of narrow absorption line s in their 
spectra (jGray &: Nagarl Il985l : ISchrijver &: Polsl Il993l : Ido Nascimento et al.l l2000l ). Several 
Doppler surveys have focused on e volved, intermediat e -mass stars on th e red giant branch 
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Hatzes et al.l 120051 ; lLovis et al.l 120051 ) and clump giant 



branch ( Sato et al. 20031 ; Setiawan et al. 2003 ). These surveys have resulted in the discovery 



of 6 substellar companions orbiting former A-type stars (Tabled]). That none of these plan- 
ets would have been detectable during their host stars' main-sequence phases highlights the 
important role evolved stars play in the study of planets around intermediate-mass stars. 

Here we present three new planet candidates around stars with M* > 1.6 M Q . These 
detections come from our precision Doppler survey of evolved stars on the s ubgiant branch 
of the H -R diagram. We discussed the selection criteria of our target stars in iJohnson et al. 
( 12006al ). along with the discovery of an eccentric hot Jupiter orbiting the 1.28 M Q subgiant 
HD 185269. We discuss our spectroscopic observations and Doppler measurement technique 
in § [21 In § |3l we present the characteristics of the host stars along with the orbital so- 
lutions for their planet candidates. We conclude with a comparison of the semimajor axis 
distributions of planets around A-type stars and lower-mass stars in § HI 



2. Observations 



We are monitorin g a sample of 159 evolved stars at Lick and Keck Observatories 
( IJohnson et al.l l2006al ). At Lick Observatory, the Shane 3 m and 0.6 m Coude Auxiliary 
Telescopes (CAT) feed the Hamilton spectrometer (jVogtlll987l ). which has a resolution of 
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Fig. 1. — Radial velocity time series for four stable subgiants with B — V colors and absolute visual 
magnitudes similar to those of the three subgiant planet host stars, a represents the standard deviation of 
the velocities about the mean (dashed lines). Observations from Keck (left panels) and Lick (right panels) 
show that subgiants in this region of the H-R diagram are typically stable to well within 10 m s _1 over time 
scales of many years. 



R ~ 50, 000 at A = 5500 A. Spectroscopic observations at Keck Observatory we re obtained 
using the HIRES spectrometer with a resolution of B, ~ 80, 000 at A = 5500 A ( jVogt et al. 



1994J ). Do ppler shifts a r e mea sured fro m each spectrum using the iodine cell method de- 



scribed by iButler et al.l (119961 ) (see also iMarcy fc Butler! 119921 ). A temperature-controlled 
Pyrex cell containing gaseous iodine is placed at the entrance slit of the spectrometer. The 
dense set of narrow molecular lines imprinted on each stellar spectrum from 5000 to 6000 A 
provides a robust wavelength scale for each observation, as well as information about the 
shape of the spectrometer's instrumental response. 
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Traditionally, the Doppler shift of each stellar observation is made with respect to an 
observed, iodine-free stellar template spectrum. These template observations require higher 
signal and resolution than normal radial velocity observations, which leads to increased 
exposure times. Given our large target list and the small aperture of the CAT, obtaining 
an observed template for each star would represent a prohibitive cost in observing time. 
We therefore perform a preliminary analysis of each star's observa tions using a synthetic , 
"morphed" template spectrum following the method described by iJohnson et al.l (j2006bl ). 
Stars showing conspicuous Doppler variations are reanalyzed using a traditional, observed 
template to verify the signal and search for a full orbital solution. 

Doppler measurements from Keck and Lick Observatories for four stable subgiants are 
shown in Figured! The error bars represent the internal uncertainties of each measurement, 
which are approximated by the weighted standard deviation of the mean velocity mea sured 



from each of the 700 individual 2 A wide chunks in each spectrum ( jButler et al.lll996l ). We 
typically achieve internal measurement uncertainties of 1-2 m s -1 for Keck observations and 
3-5 m s -1 at Lick. Subgiants have an additional 4-6 m s _1 of "jitter" — velocity scatter 
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value of 5 m s 1 for our subgiants, which is added in quadrature to the internal uncertainties 
of the measurements before searching for a best-fit orbital solution. 

After determining the best-fit Keplerian solution using a Levenberg-Marquardt, least- 
squares minimization, we estimate the orbital parameter uncertainties using a bootstrap 
Monte Carlo method. We first subtract the best-fit Keplerian from the measured velocities. 
The residuals are then scrambled and added back to the original measurements, and a new set 
of orbital parameters is obtained. This process is repeated for 1000 trials, and the standard 
deviations of the parameters from all trials are adopted as the formal, la uncertainties. 



3. Stellar Properties and Orbit Solutions 

3.1. Estimates of Stellar Properties 

We estimated the stellar properties of our target stars using two pr i mary methods: 
the LTE spectral synthesis method (SME) described by IValenti fc Fischer! (120051 ). and the 
Padovei stellar interior models. The spectral synthesis method uses a non-linear least- 
squares algorithm to vary the parameters of a synthetic spectrum to search for a fit to an 



2 See also 



http: / /pleiadi. pd.astro.it / 
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Fig. 2. — This figure illustrates the interpolation method employed to determine accurate stellar masses for 
the planet host stars. We estimated each star's mass, radius and age by interpolating its Hipparcos B—V color 
and absolute visual magnitude My onto grids of four different metallicities: [Fe/H] = —0.4, 0.0, +0.18, +0.3 
(open diamonds). For each star's measured value of B — V and My, the three lines show the dependence 
of stellar mass on the measured [Fe/H], estimated using a cubic spline interpolation between the diamonds. 
Similar dependencies were determined for stellar radii, luminosities and ages. 



iodine-free stellar template spectrum. The free parameters in the fit are the abundances 
of heavy elements; effective surface temperature, T eff ; surface gravity, log q; and broad ening; 
effects due to the star's projected rotation velocity, V^otsinz. IValenti fc Fischer! (120051 ) esti- 
mate a precision of 0.04 dex in metallicity, 44 K in effective temperature, 0.3 dex in logg, 
and 0.5 km s _1 in rotational velocity. 

To estimate stellar masses, radii, luminosities and ages, we used the Padova theo- 
retical stellar mode ls, which have been transformed into several photometric systems by 
Girardi et al.l (120021 ). Stellar properties can be inferred by interpolatin g a star's color , abso- 
lute magnitude and metallicity onto these model grids. However, the iGirardi et al.l (120021 ) 
model grids are defined at widely-spaced metallicity intervals, with [Fe/H] = -0.4, 0.0, +0.18 
and +0.30. Since the uncertainties in our spectroscopically derived metallicity estimates are 
much less than the model grid intervals, and because the derived stellar properties do not 
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vary linearly with [Fe/H], we could not simply perform a linear, 3- dimensional interpolation 
of My, B — V and [Fe/H]. Instead, we first linearly interpolate the stars' colors and absolute 
magnitudes onto each of the four metallicity grids. We then use a cubic spline interpolation 
between the grid points to measure the desired stellar property (e.g. mass) at the star's 
measured [Fe/H]. Our procedure is illustrated in Figure [21 which shows stellar mass as a 
function of [Fe/H] for each star's absolute magnitude and color. The same procedure was 
used for stellar radii, luminosities and ages. 



We compared our interpolated stellar properties to the iTakeda et al.l (120071 ) theoretical 
interior models of the s tars in the Spectroscopic Properties of Cool Stars catalog (SPOCS 
Valenti &: Fischer! 120051 ). We found a subset of 11 evolved stars in the catalog with 2.0 < 
My < 3.0 and 0.7 < B — V < 1.1. Differences between our inferred values and those 



from ITakeda et al.l (120071 ) had an rms scatter of 7% in mass, 12% in radius, with a median 
offset of -2% and -4% in each parameter, respectively. Ages of this subset of evolved stars 
estimated by the two methods have a difference of -0.4 Gyr with and rms scatter of 1.1 Gyr. 
We therefore adopt fractional uncertainties of 7% for our derived masses, 12% for radii and 
1 Gyr for ages. We list the full set of derived stellar properties of the three candidate 
planet host stars in Table [2J We summarize each star's properties and orbital solution in 
the following subsections. 



3.2. HD 192699 



HD 192699 (HIP 99894) is listed with a G5 spectral type in the Hipparcos Catalog, with 
V = 6.44, B — V = 0.867 and a parallax-based distance of 67 pc (jESAlll997l ). However, no 
luminosity class is given. Based on its distance, we calculate My = 2.30, which at its B — V 
color places the star 3. 7 mag a bove the mean main-sequence of stars in the Solar neigh- 
borhood, as defined by IWrightl (12004 ). Based on its color and absolute magnitude, we find 
that HD 192699 is likely a G8 IV subgiant near the base of the red giant branch. Commen- 
surate with its evolved status, HD 192699 is chromospherically inactive, with S = 0.12 and 



HK 



5.29 as me asured from the Call H&K line core and averaged over all observations 



(I Wright et al.ll2004h . 



Based on our LTE spectral analysis, we find that HD 192699 is metal-poor, with [Fe/H] = 
—0.15, and slowly rotating, with V rot sm.i = 1.9 km s _1 . The other stellar parameters derived 
from our spectral analysis are liste d in Table El We in terpolated the star's color, absolute 
magnitude and metallicity onto the lGirardi et al.l (120021 ) theoretical stellar model grids using 
the method described in § 13.11 Our interpolation yields a stellar mass M* = 1.68 M , radius 
R* = 4.25 R & , and an age of 1.8 Gyr. 
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Fig. 3. — Radial velocity time series for HD 192699 measured at Lick Observatory. The error bars reflect 
the quadrature sum of the internal measurement uncertainties and 5 m s _1 of jitter. The dashed line shows 
the best-fit orbital solution, which has \[y^= 1.14. 



We began observing HD 192699 in 2004 May at Lick Observatory using the 3 m Shane 
Telescope and 0.6 m CAT. Table [3] lists our 34 velocity measurements, along with their 
times of observation and internal measurement uncertainties (without jitte r). Our first 7 



observations, initially analyzed using a synthetic stellar template spectrum (IJohnson et al. 



2006bl ). showed correlated variations spanning two observing seasons. We obtained a high- 
quality observed template using the Shane 3m telescope and initiated intensive follow-up 
observations during the Fall 2006 observing season. The Keplerian signal is visible to the 
eye (Figure [3]), obviating a periodogram analysis. 

The best-fit Keplerian orbit has a period of P = 351.5 d, velocity amplitude K = 
51.3 m s _1 , and eccentricity e = 0.149 ± 0.06. With an assumed stellar mass of 1.68 M , we 
estimate a minimum planet mass Mp sin i = 2.5 Mj up and orbital separation a = 1.16 AU. 
The fit has rms = 9.2 m s _1 and a reduced \fxl = 1.14, consistent with the measurement 
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errors and jitter. The full set of orbital parameters and uncertainties is listed in Table HI 



3.3. HD 210702 



HD 210702 (HIP 109577, HR8461) is listed in the Hipparcos catalog as a Kl star (no 
luminosity class given) with V = 5.93, B — V = 0.951, and a parallax-based distance of 56 pc. 
Given its distance and apparent magnitude, we calculate an absolute magnitude My = 2.19, 
which places i t 4.2 mag above the average main sequence of stars in the Solar neighborhood 
( iWrightl I2004J ) . We therefore estimate that HD 210702 is a class K1IV subgiant near the 
base of the red giant branch. 

Based on our LTE spectral analysis, we find that HD 210702 is somewhat metal-rich, 
with [Fe/H] =+0.12±0.04, and slowly rotating, with V^. ot sini = 1.9 k m s~ l . Our interpo- 
lation of the star's color, absolute magnitude and metallicity onto the iGirardi et al.l (120021 ) 
stellar model grids yields a stellar mass = 1.85 M & , stellar radius = 4.72 R Q , and an 
age of 1.4 Gyr. Consistent with its post-main-sequence evolutionary status, HD 210702 is 
chromosp herically inactive w ith S = 0.11 and R'hk = -5.35, as measured from its Call H&K 
emission ( IWright et al.ll2004l ). The other stellar parameters derived from our spectral anal- 
ysis and stellar model interpolation are listed in Table [2J 

We began monitoring HD 210702 in 2004 August at Lick Observatory. The first 9 obser- 
vations were Doppler-analyzed using a synthetic template, and showed excessive variability 
with rms = 19 m s _1 . We then obtained a traditional, observed template to confirm the vari- 
ations with higher Doppler precision. The full set of velocities is listed in Table [6] (without 
jitter) and plotted in Figure HI The error bars in Figure H] have been augmented by adding 
5 m s _1 of jitter in quadrature to the internal measurement uncertainties. 

The best-fit Keplerian orbital solution is shown in Figure HJoverplotted on the velocities. 
The solution has a 341.1 day period, an eccentricity e = 0.152 ± 0.08, and a semiamplitude 
K = 39.2 m s _1 . The fit residuals have rms = 7.4 m s _1 and reduced v/xE = 1-07, consistent 
with the internal measurement uncertainties and jitter. Assuming a stellar mass M* = 
1.85 Mq, the best-fit solution yields a relative separation a = 1.17 AU. 

We find that the inclusion of a linear trend in the orbital solution yields a slight im- 
provement in the quality of fit, decreasing the rms scatter of the residuals from 7.4 m s _1 
to 6.7 m s -1 , and the reduced \fxl from 1.07 to 1.00 after accounting for the extra free 
parameter in t he Keplerian- p lus-tr end model. We tested the validity of the trend using the 
prescription of IWright et al.l (120071 ) . and found a false-alarm probability of 49%. The large 
FAP indicates that the apparent linear trend is likely due to noise rather than and additional 



CO 

e 



o 

1) 
> 

"3 



50 - 

; 

o -! 



& -501- 



o 



25 F 


-25L 



-10- 



MpSini = 2.0 M Jup 
P = 341.1 d 
K = 39.2 m s" 1 
e = 0.15 ^ 



HD 210702 
Lick Obs. 



* 



L' x 



f 



+ 



s / 



-1 



a = 7 A m s 
I ■ ■ ■ I ■ ■ ■ I 



v 



T 



200 400 600 800 1000 
JD - 2453200 

Fig. 4. — Radial velocity time series for HD 210702 measured at Lick Observatory. The dashed line shows 
the best-fit orbital solution, which has \f\i= 1.07. 



orbital companion. Indeed, the trend appears to be driven primarily by the three outliers 
near JD = 100, 400 and 800 (Figure HJ). We therefore favor the single-planet Keplerian 
model summarized in Table HJ 



3.4. HD 175541 



HD 175541 (HIP 92895) is listed in the Hipparcos Catalog asa G8V star with V = 8.03, 
B — V = 0.869 and a parallax-based distance of 128 pc (jESAlll997l ). Given its distance, the 
star has My = 2.49, placing it 3.5 mag above the mean main-sequence of stars in the Solar 
neighborhood (jWrightll2004l ). Like most evolved stars, HD 175541 is chromospherically quiet 
with S = 0.11 and R' HK = -5.28 (jWright et al.ll2004h . Its low chromospheric activity and 
location in the H-R diagram indicate that HD 175541 is most likely a luminosity class IV 



GO 



40 



20 - 



- 



o 
o 

> 

1 -201- 

(0 

02 



-40 



HD 175541 
Keck Obs. 







(7 = 11.6ms 



f ♦ 



-i 



1 000 2000 3000 

JD - 2450300 



Fig. 5. — Radial velocities for HD 175541 measured at Keck Observatory. The error bars represent the 
internal uncertainty of each measurement without accounting for stellar jitter. 



star on the subgiant branch, rather than a class V dwarf. 



HD 175541 is listed in the SPOCS Catalog (IValenti fc Fischer! 120051 ) with a metal abun- 
dance slightly below Solar ([Fe/H] =-0.07±0.04) and projected rotational velocity V rot smi = 
2.9 km s " 1 . Interpola t ion of the star's B — V color, absolute magnitude and metallicity 
onto the iGirardi et al.l (120021 ) stellar model grids yields a stellar mass = 1.65 ra - 
dius i?* = 3.85 R Q , and an age of 1.9 Gyr. The interior models of iTakeda et al.l (120071 ) 
yield M* = 1.5 2 M^and R* = 3 72 E^ . The SPOCS Catalog lists M* = 1.74 M , and 
= 4.11 Rq (jValenti fc Fischer! 120051 ). The variances of these different mass and radius 
estimates are 0.1 M & and 0.2 R Q , respectively, which are consistent with our estimate of 
uncertainties in §3.11 The other stellar properties are listed in Table [2j 



-12- 



121 — ' — 




10 100 1000 

Period [days] 




Fig. 6. — Top: Periodogram analysis of the RV time series of HD 175541. A strong peak is visible near 
P = 300 d with an analytic false-alarm probability FAP< 0.1%. Bottom: Empirical assessment of the FAP of 
the best-fit Keplerian model. The original, unscrambled velocities yield an orbital solution with \[~X%,= 1.01 
(arrow). The histogram shows the distribution of yxj obtained from the best-fit orbital solution for each 
of the scrambled-velocities trials. None of the 10 4 trials produced a value of •v/xj lower than the value 
obtained from the original time series, resulting in FAP< 0.01% (cf § I3.4|) . 
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Fig. 7. — Radial velocity observations of HD 175541 phased at P = 297.3 d. The gray points lie outside of 
phases 0.0 and 1.0 and are included to guide the eye. The dashed line shows the best-fit Keplerian orbital 
solution, which has \/xf = 1.01. 



HD 175541 was one of the original stars added to the CCPS Keck program in 1996, 
and was subsequently added to our list of intermediate-mass stars in 2004. Table [5] lists 
our 29 Doppler measurements along with their observation dates and internal uncertainties 
(without jitter). Figure shows that the rms scatter of the velocity measurements is a factor 
of 6 greater than the mean internal uncertainty (<f v fa 2 m s _1 ), and 2-3 times larger than 
the rms scatter of stable Keck subgiants (Figured]). A Lombe-Scargle periodogram analysis 
of the velocities reveals a pronounced peak near P = 300 d, with an analytical false-alarm 
probability < 0.1% (Figure E]). 

To search for the best-fit orbital solution, we added 5 m s _1 of jitter in quadrature 
to the internal measurement uncertainties. We find that a Keplerian with P = 297.3 d, 
K = 14.0 m s -1 and e = 0.33 provides the best fit to the data, resulting in rms = 5.6 m s~ x 




= 1.01. Figure H shows the the radial velocities phased at P = 297.3 day, along 
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with the best -fit orbital solution (the gray points show the measurements at phases outside 
of phases 0.0 and 1.0, in order to guide the eye). Assuming a stellar mass of 1.65 M , we 
estimate a minimum planet mass Mpsinz = 0.61 Mj up and orbital separation a = 1.03 AU. 

While the strong periodogram peak and low v/xf are indicative of a correlated signal 
resulting from an orbiting planet, it is still possible that random variability could conspire to 
produce a false periodicity in our sparse series of measurements. To test the nul l hypo thesis, 
we used the "scrambled" velocity false-alarm test described by iMarcy et al.l (120051 ). For 
10 4 separate trials, we held the observation times constant and scrambled the order of the 
measurements using a pseudo random number generator. This has the effect of keeping 
the sampling constant while removing any true temporal coherence, if such a signal exists. 
For each of the scrambled trials, we perform a full search for the best-fit Keplerian orbital 
solution — with jitter — and record the \[y? v from the fit. 

The distribution of v/xf generated from the scrambled-velocity trials is then compared 
to the fit obtained from the original time series, as shown in the lower panel of Figure [6j 
None of the 10 4 scrambled trials produced a a/xT equal to or lower than the best-fit solution 
to the original time series, resulting in a false-alarm probability of < 0.01%. From this 
test, we conclude that the temporally correlated signal seen in the velocity time series is 
likely real, rather than an artifact of random noise. We find that the best explanation of the 
periodic signal is the presence of an unseen planetary companion orbiting HD 175541. 



4. Summary and Discussion 

We present precision Doppler measurements of three intermediate-mass subgiants that 
show periodic variations in their radial velocities consistent with planet-mass orbital com- 
panions. Inte r polat ion of the stars' absolute magnitudes, colors and metallicities onto the 



Girardi et al.l (120021 ) stellar interior models shows that all three stars have masses ranging 
from 1.65 M & to 1.85 M & . Figure [H] shows these massive host stars on an H-R diagram, 
along with their theoretical evolution tracks. Following the tracks back to the zero-age main 
sequence reveals that these present-day subgiants were originally early-type dwarfs with 
B — V < 0.2 and spectral types ranging from A2V to A5V. The three long-period planets 
presented here would not have been detectable during their stars' main-sequence phases 
due to the jitter and rotational line broadening typical for intermediate-mass dwarfs. These 
planets orbiting "retired" A stars illustrate how evolved stars provide a unique window into 
stellar mass and planetary domains otherwise inaccessible to Doppler-based planet searches. 

There are now 9 former A-type stars (1.6 < M* < 3.0 M ) with planetary companions. 




Fig. 8. — H— R diagram illustrating the properties of the three subgiant planet host stars (pentagrams) 
compared to their main-sequence progenitors (filled circles). The connecting lines represent each star's 
Girardi et all \200j \ theoretical mass track, interpolated for that star's metallicity. The thick, diagonal line 



is the theoretical zero-age main sequence for [Fe/H]=0.0. The hashed region shows the approximate range 
of colors and magnitudes of stars with spectral types A2V-A8V. 



We list some of the properties of these massive host stars and their planets in Table [TJ All 
9 planets orbit beyond ~0.78 AU from their stars. This paucity of planets with semimajor 
axes a < 0.78 AU is unlikely to be due to a detection bias. For a given planet mass and 
stellar mass, the velocity semiamplitude of a star scales as K ~ a" 1 / 2 , making planets in 
smaller orbits easier to detect. The detectability of close-in planets is also facilitated by the 
increased number of orbital cycles that are observable over a given time span. 

We consider two possible explanations for the observed lack of close-in planets around 
intermediate-mass stars. The first possibility is that planets around A-type stars have 
the same semimajor axis distribution as planets orbiting lower-mass stars, but the close-in 
planets were destroyed by the expanding atmospheres of their giant host stars. Alternatively, 
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planets orbiting A-type stars may have a different semimajor axis distribution than lower- 
mass stars, with planets residing preferentially in long-period orbits beyond ~ 0.8 AU. 

These possibilities can be explored by comparing the properties of planets in Tabled] to 
planets orbiting lower-mass stars listed in the CNE. We exclude extremely low-mass planets 
with K < 15 m s _1 that would not be easily detectible around higher-mass subgiants and 
giants. We use a one-sided Kolmogorov-Smirnov (K-S) test to compare the semimajor 



axis distributions of planets around intermediate-mass and lower-mass stars (IPress et al. 



19921 ). We find the probability that the two distributions are identical is only 0.06%. Under 
the assumption that the semimajor axis distribution of planets is independent of stellar 
mass, short-period planets orbiting evolved A-type stars must be efficiently destroyed by 
the expanding atmospheres of their giant host stars. The validity of this hypothesis depends 
on whether the radii of ~ 2 M Q giants are large enough to engulf planets out to ~ 0.8 AU. 



Figure M shows the evolution of the radius of a 2.0 M & star according to the lGirardi et al. 



(120021 ) stellar evolution models. As the star crosses the Hertzsprung Gap during its subgiant 
phase, its radius remains nearly constant at a = 5 R Q = 0.023 AU, which is within the orbit 
of a P = 3 day hot Jupiter. Not until the star begins to ascend the RGB does its outer 
atmosphere begin to encroach on the orbits of short-period planets. But even at the tip of 
the RGB (near the helium flash), the radius of a 2 M Q star is only at the distance of a 10 day 
hot Jupiter at a ~ 26 R & =0.12 AU (the radius of a 2.5 M & red giant is not much larger at 
a ?s 32 Rq = 0.15 AU). Thus, engulfment cannot be solely responsible for the lack of close-in 
planets around subgiants and K giants. Indeed, engulfment can only be important for 4 of 
the stars in Table [TJ the post-helium-flash clump giants HD 104985, HD 11977 and e Tau, 
and HD 13189 which has a poorly constrained radius due to its highly uncertain parallax. 

The evolution of planetary orbits from 0.05-0.15 AU in the presense of an expanding 
stellar atmosphere has not been examined in d etail. The effects of planet engulfment on its 



host star have been studied by ISiess fc Livid (119991 ) . but a key assumption in their model 



is that the substellar companion is destroyed. Since it is unclear what happens to a planet 
when it interacts with the atmosphere of its expanding host star, we simply assume that 
planets orbiting within the radius of a giant star are destroyed^]. Under this assumption, 
we would expect a deficiency of hot Jupiters around clump giants out to ~0.15 AU, but no 
corresponding deficiency around subgiants and K giants. 

We now analyze the lack of close-in planets around the sample in Table [1] accounting 



Maxted et all |2006t) discovered a short-period substellar companion that apparently survived engulf- 



ment as its parent star evolved into a white dwarf. However, no Jovian planet has yet been detected around 
a white dwarf. 
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Fig. 9. — Radius of a 1.9 M & star as it evolves off of the main sequence (becoming redder). 
The horizontal dashed line s depict the semimajor axes of planets with periods of 3 and 10 
days. The vertical dot-dashed line shows the approximate B — V color of the star as it 
begins to ascend the red giant branch. Because the radii of subgiants are small enough to 
avoid interference with close-in planets, our Doppler survey is sensitive to the same range 
of orbital separations as surveys of main-sequence stars. Planets with a < 30 R Q may be 
destroyed by the expansion of their host stars on the red giant branch. 
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the possible destruction of hot Jupiters around clump giants. For subgiants and giants 
we can use the K-S test as before, which yields a probability of 0.7% that the semimajor 
axis distribution is the same as lower-mass stars in the CNE. For clump giants we exclude 
planets from the CNE with a < 0.15 AU, and the corresponding probability from the K-S 
test is 1.7%. Thus, the distribution of close-in planets around former A-type stars remains 
inconsistent with the distribution of planets in the CNE. Since engulfment does not provide 
an adequate explanation for the lack of close-in planets in Table [TJ we are left with the 
possibility that the semimajor axis distribution of planet around A-type stars is significantly 
different than the distribution around lower-mass stars (M* < 1.6 M Q ). 

Differences between the se mimajor axes of plan ets around stars of various masses has 



previously been investigated by lBurkert & Idal (120061 ). From their study of the orbital prop- 



erties of known exoplanets, they find evidence of a gap in the semimajor axis distribution 
around stars with masses M* > 1.2 M , with fewer planets between 0.08 AU and 0.6 AU 
compared to lower-mass stars. They were able to reproduce this gap in their Monte Carlo 
simulations of planet migration, and they attribute the gap to the shorter depletion timescales 
of disks around intermediate-mass stars. 

The semimajor axis distribution of planets as a function of stellar mass can be in- 
vestigated further with the inclusion of a larger sample of intermediate-mass subgiants in 
Doppler-based planet searches. As Figure [9] shows, Doppler surveys of subgiants can probe 
occurrence of Jovian planets at orbital distances ranging from many AU down to as close 
as 0.05 AU, the realm of hot Jupiters. The smaller radii of subgiants also result in higher 
surface gravit i es com pared to giants, which leads to lower levels of pulsation-induced jitter. 



Hekker et al.l (120061 ) show that giants with B — V > 1.2 typically have jitter values greater 
than 20 m s _1 , ostensibly due to radial and non-radial pulsation modes. Only giants blue- 
ward of this limit are stable to within 20 m s" 1 , compared to the 4-6 m s" 1 of jitter seen in 
subgiants (B — V< 1.0, My < 2.0). This increased velocity stability, coupled with their rela- 
tively small radii, therefore make subgiants ideal proxies for A-type dwarfs in Doppler-based 
planet searches. 

The primary limitation of subgiants is their relative scarcity, which restricts the num- 
ber of bright targets suitable for high-resolution spectroscopic observations. The time it 
takes stars to cross the Hertzsprung Gap is small compared to the star's lifetime — of order 
100 Myr — rendering Hertzsprung Gap stars within 200 pc rare compared to main-sequence 
stars and giants. Additional targets can be found further from the Sun, with fainter apparent 
magnitudes (V > 7.5). In the near future, we plan to expand our sample of subgiants using 
the Keck telescope and HIRES spectrometer in order to further investigate the orbital prop- 
erties, planet masses and occurrence rate of planets orbiting intermediate-mass stars. As the 
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number of subgiants included in Doppler surveys increases, it will become apparent whether 
the lack of short-period planets around intermediate-mass stars is a result of different for- 
mation and migration mechanisms in the disks of A-type stars, or simply a consequence of 
the small number of massive subgiants currently surveyed. 
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Table 1. Planet Host Stars With M* > 1.6 M a 



HD Spectral Type M* (M ) i?* (R e ) a (AU) References 



13189 


K2 II 


2-6 b 




1.5-2.2 


1 


28305 c 


KOIII 


2.7 


13.7 


1.93 


2 


11977 


G5III 


1.9 


13 


1.93 


3 


62509 d 


KOIII 


1.86 


8.8 


1.7 


4,5 


210702 


K1IV 


1.85 


4.72 


1.17 


6 


192699 


G8IV 


1.68 


4.25 


1.16 


6 


175541 


G8IV 


1.65 


3.85 


1.03 


6 


89744 


F7IV 


1.65 


1.1 


0.93 


7 


104985 


G9III 


1.6 


8.9 


0.78 


8 



Excluded from this list of evolved, intermediate-m ass planet host stars 



are 7 Cep A and HD 47536. The discove ry paper by 



jatzes et all (120031 ) 



cites a stellar mass of 1.59 M & . However, iTorred (120071 ) find a much lower 
value of M* = 1. 181^11 -^0 by usin g a theoretical model gri d interpolation, 



and a recent dynamical analysis by lNeuhauser et all (120071 ) yields a stellar 
mass of 1.40 ± 0.12 M . For the clump giant HD 47536, ISetiawan et al 



(120031 ) list a highly uncertain mass with a lower limit of 1.1 M , which 
falls well below the 1.6 M cutoff used for this table. Also omitted are the 



3 intermediate-mass K giants in iMitchelll (120041 ). which have not yet been 
published in a refereed journal. 



Schuler et all (120051 ) 



= HD 28305 = e Tau. 



d HD 62509 = Gem = Pollux. 



References. — (1 ) lHatzes et all (I2005h: [ 21 ISato et all (1200711: (3 ) 
Setiawan et aD (l2005h: (A) lHatzes et al.l2006h: fsTlReffert et al.l J2006h: 
(6) this work; (7) korzennik et all J2OO0I ): (8) ISato et al.l j200ah ~ 
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Table 2. Stellar Parameters 



Parameter 


HD 192699 


HD 210702 


HD 175541 


V 


6.44 


5.93 


8.03 


My 


2.30 


2.19 


2.49 


B-V 


0.867 


0.951 


0.869 


Distance (pc) 


67 


56 


128 


[Fe/H] 


-0.15 (0.04) 


+0.12 (0.04) 


-0.07 (0.04) 


r e// (K) 


5220 (44) 


5010 (44) 


5060 (44) 


K-ot sini (km s _1 ) 


1.9 (0.5) 


1.7 (0.5) 


2.9 (0.5) 


logs' 


3.44 (0.3) 


3.29 (0.3) 


3.52 (0.3) 


M, (M ) 


1.68 (0.12) 


1.85 (0.13) 


1.65 (0.12) 


-R* (-^©) 


4.25 (0.51) 


4.72 (0.57) 


3.85 (0.46) 


I/* (£©) 


11.5 (0.3) 


13.1 (0.3) 


9.56 (0.3) 


Age (Gyr) 


1.8 (1.0) 


1.4 (1.0) 


1.9 (1.0) 




0.12 


0.11 


0.11 




-5.29 


-5.35 


-5.28 
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Table 3. Radial Velocities for HD 192699 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



13155. 


.988 


-60.66 


5.14 


13257. 


.741 


3.44 


6.31 


13522. 


.921 


-47.84 


4.90 


13576. 


901 


-20.91 


5.33 


13619. 


.774 


13.49 


4.90 


13640. 


.713 


26.10 


4.99 


13668, 


.641 


19.02 


4.66 


13879. 


.905 


-46.41 


4.39 


13891. 


.854 


-34.19 


5.46 


13894. 


.944 


-42.03 


4.94 


13895. 


.904 


-56.83 


4.71 


13896. 


937 


-50.29 


4.57 


13905. 


962 


-28.97 


6.96 


13921. 


.932 


-18.30 


5.37 


13922. 


.851 


-1.89 


4.63 


13926. 


.856 


-18.51 


4.47 


13951. 


.801 


11.70 


4.62 


13959. 


.771 


-4.84 


5.15 


13966. 


.784 


8.99 


5.63 


13975. 


.738 


5.99 


4.65 


13976. 


.714 


14.68 


4.87 


13998. 


.725 


25.11 


5.28 


14001. 


.749 


34.99 


4.68 


14020. 


.654 


49.36 


5.02 


14021. 


.682 


25.41 


5.42 


14034. 


.638 


18.33 


5.23 


14035. 


.671 


9.67 


5.46 


14039. 


603 


-2.60 


5.69 


14046. 


.658 


6.44 


5.05 


14049. 


615 


5.84 


5.60 
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Table 3 — Continued 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



14059.610 -5.29 4.55 

14070.664 -16.58 5.30 

14072.612 -18.70 5.17 

14170.058 -78.23 5.92 
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Table 4. Orbital Parameters 



Parameter 


HD 192699 b 


HD 210702 b 


HD 175541b 


P(d) 


351.5 (6) 


341.1 (7) 


297.3 (6) 


V (JD) 


2452994 (30) 


2453118 (40) 


2450213 (20) 


e 


0.149 (0.06) 


0.152 (0.08) 


0.33 (0.2) 


Ki (m s- 1 ) 


51.3 (5) 


39.2 (4) 


14.0 (2) 


lu (deg) 


54 (30) 


301 (30) 


183 (30) 


M P sin % (M Jup ) 


2.5 


2.0 


0.61 


a (AU) 


1.16 


1.17 


1.03 


Fit RMS (m s- 1 ) 


9.2 


7.4 


5.6 


-\ / 
V ^-v 


1.14 


1.07 


1.01 


N obs 


34 


29 


29 



a Time of periastron passage. 
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Table 5. Radial Velocities for HD 175541 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



10283. 


.952 


5.14 


1.46 


10603. 


.027 


10.95 


1.63 


10665. 


908 


12.81 


1.63 


11011. 


.878 


0.94 


2.19 


11069. 


.830 


-14.86 


2.17 


11367. 


.894 


-7.68 


1.98 


11410. 


.857 


-12.21 


2.29 


11441. 


.759 


-5.03 


1.96 


11705. 


.979 


-15.48 


2.25 


11984. 


.164 


-13.59 


2.50 


12004. 


.126 


-15.69 


2.18 


12008. 


065 


-19.07 


2.17 


12009. 


.118 


-28.52 


2.44 


12030. 


989 


-4.26 


1.67 


12095. 


962 


5.95 


1.94 


12098 


016 


-2.49 


2 09 


12099. 


039 


10.75 


1.71 


12100. 


.956 


2.62 


1.94 


12127. 


.874 


5.87 


2.05 


12391. 


.132 


3.80 


2.13 


12488. 


.898 


10.77 


1.78 


12573. 


.725 


-9.38 


1.82 


12833. 


.948 


-1.19 


1.82 


13239. 


.787 


0.00 


1.80 


13546. 


.914 


7.36 


1.98 


13603. 


.852 


15.72 


1.89 


13807. 


.150 


-3.79 


1.93 


13932. 


.962 


16.87 


1.76 


13968. 


.920 


12.16 


1.83 
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Table 6. Radial Velocities for HD 210702 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



13241. 


.863 


18.13 


3.32 


13256. 


.776 


9.30 


3.91 


13310. 


.752 


-16.07 


4.33 


13342. 


.598 


-0.18 


3.37 


13569. 


930 


42.99 


4.77 


13573. 


.904 


29.79 


5.13 


13574. 


.927 


38.71 


4.47 


13575. 


.897 


29.80 


4.41 


13640. 


.779 


10.99 


5.45 


13922. 


930 


28.01 


3.39 


13959. 


.860 


0.18 


3.22 


13975. 


.798 


-1.47 


3.54 


13976. 


.764 


-8.51 


3.20 


13998. 


.768 


-10.73 


3.52 


14001. 


.789 


-13.09 


3.61 


14020. 


721 


-15 97 

la. j 1 


3 73 


14021. 


.735 


-11.53 


4.38 


14034. 


.695 


-22.55 


4.41 


14046. 


.694 


-31.56 


3.40 


14059. 


673 


-14.62 


3.57 


14070. 


.705 


-11.10 


3.70 


14092. 


.594 


1.07 


4.06 


14097. 


.649 


-2.40 


3.31 


14098. 


603 


1.90 


3.78 


14103. 


.598 


9.35 


3.65 


14109. 


.626 


9.02 


5.12 


14117. 


.595 


22.46 


5.58 


14136. 


600 


38.79 


3.09 


14197. 


,033 


53.95 


4.26 
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